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ABSTRACT 

Water (H 2 0) ice is an important solid constituent of many astrophysical environments. 
To comprehend the role of such ices in the chemistry and evolution of dense molecular 
clouds and comets, it is necessary to understand the freeze-out, potential surface reac- 
tivity, and desorption mechanisms of such molecular systems. Consequently, there is 
a real need from within the astronomical modelling community for accurate empirical 
molecular data pertaining to these processes. Here we give the first results of a labo- 
ratory programme to provide such data. Measurements of the thermal desorption of 
H2O ice, under interstellar conditions, are presented. For ice deposited under condi- 
tions that realistically mimic those in a dense molecular cloud, the thermal desorption 
of thin films («50 molecular layers) is found to occur with zero order kinetics charac- 
terised by a surface binding energy, Ed es , of 5773 ±60 K, and a pre-exponential factor, 
A, of 10 30±2 molecules cm~ 2 s -1 . These results imply that, in the dense interstellar 
medium, thermal desorption of H2O ice will occur at significantly higher temperatures 
than has previously been assumed. 
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1 INTRODUCTION 

It has become clear over the last decade that purely gas 
phase schemes cannot account for the variety and rich- 
ness of chemistry occurring in the interstellar medium 
(ISM), particularly in denser regions of molecular clouds 
where star formation is occurring. In such environ- 
ments gas-grain interactions must also play a key role 
(Williams 1998). Within the lifetime of a dense molec- 
ular cloud (on the order of 10 6 yr), the grains accrete 
icy mantles as atomic and molecular species freeze-out 
onto the grains. The ices are processed further by cosmic 
ray impacts, UV photolysis, and shocks. Consequently, 
the chemical composition of these icy mantles is signif- 
icantly different from that of the local gas phase en- 
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vironment (Rawlings 2000). During cloud collapse and 
stellar formation the grains are reheated, and molecules 
are recycled from the solid state into the gas phase envi- 
ronment. In regions where the H2O concentrations are 
depleted it is assumed that the species is still locked in 
the icy mantles on the grains, and consequently, that 
the temperature is below the sublimation temperature 
of the tracer gas (van Dishoeck & van der Tak 2000). 
When evaluating these environments it is clearly neces- 
sary to understand the sublimation behaviour of the ice, 
and its thermal stability. The surface binding energy, 
Erf es , and sublimation rates of molecular ices may only 
be obtained empirically. When combined, these data 
can be used to determine molecular residence times on 
the ice surface, in terms of a population half-life as a 
function of temperature. These data are also applicable 
in other astronomical environments, such as cometary 
nuclei and comae, planetary surfaces and satellites. 

H2O is the most abundant molecular ice in dense 
ISM regions (Ehrenfreund & Schutte 2000). The H 2 
ice band at 3.07 /im is detected on lines of sight towards 
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reddened stars that have a visual extinction above a 
critical value (Whittet 1993) and is typically one of the 
strongest bands in interstellar IR spectra. The abun- 
dance and morphology of H2O ice in the ISM is gen- 
erally inferred from the intensity and line profile of the 
3.07 /im band. At low pressures and temperatures, such 
as those in the ISM, the H2O ice can exist in both 
amorphous and crystalline forms. To further compli- 
cate the picture, the amorphous ice can exist in both 
high-density (1.1 g cm~ 3 ) and low-density (0.94 g cm~ 3 ) 
forms, reflecting differences in the porosity of the mate- 
rial, with a phase change from former to latter occurring 
irreversibly between 38 - 80 K (Jcnniskcns et al., 1995). 
It is generally assumed that the dominant morphol- 
ogy of icy mantles on interstellar grains resembles that 
of high density amorphous ice, with other 'impurity' 
molecules trapped in the H2O ice matrix (Jcnniskcns et 
al. 1995). This is evident in comparisons of ISO SWS 
and ground based observations with laboratory spectra 
of mixed ice analogues (Ehrenfreund & Schutte 2000 
and references therein). However, laboratory spectra of 
pure H2O ices show clear differences in the line profile 
of the 3.07 /im feature associated with different ice mor- 
phologies and temperatures (Jenniskens et al. 1997). It 
is therefore likely that some H2O ices have undergone 
'mild processing' and re-cooling, (such as gentle heat- 
ing below 80 - 100 K, or rapid localised heating due 
to shocks), and will therefore also exist in low-density 
amorphous and crystalline forms (Jenniskens & Blake, 
1994, 1996). The physical properties of the H2O ice, 
such as density, conductivity, vapour-pressure, and sub- 
limation rate, are dictated by its structure. Significant 
differences are expected between the surface chemistry 
and bulk behaviour of the ice phases. In addition, the 
physical and chemical properties of the ice may also be 
affected by way the ice film is deposited, its lifetime, and 
processing prior to thermal desorption (Sack & Baragi- 
ola 1993). 

This paper is the first in a series of papers report- 
ing new experimental results on gas-surface interactions 
occurring under conditions resembling those in the ISM 
and star- forming regions. With the aid of a purpose- 
built instrument, we are now able, for the first time, to 
provide the astrochemical community with experimen- 
tal results addressing not only desorption energies, but 
also sticking probabilities and desorption mechanisms. 
Here we present the results of laboratory studies to de- 
termine the surface binding energies on, and sublima- 
tion rates of H2O ice under interstellar conditions. Pre- 
vious studies, where the surface binding energies were 
determined from spectroscopic data, suggested that the 
surface binding energy of H2O on H2O ice was 4815 
±15 K or 5070 ±50 K, for unannealed and annealed ice 
samples respectively (Sandford & Allamandola 1988). 
It was also presumed that the desorption kinetics of the 
system were first order. These data have been used ex- 



tensively in the astronomy literature, where it is now 
widely reported that under interstellar conditions H2O 
can only exist in the solid phase below about 100 K 
(Hasegawa & Herbst 1993, Sandford & Allamandola 
1993, Ehrenfreund & Schutte 2000). The direct tech- 
niques used to determine the surface binding energy in 
this work are described in §2. The results are detailed in 
§3. In these experiments the thermal desorption of thin 
H2O ice films is found to occur with zero order kinet- 
ics characterised by a surface binding energy, Ed es , of 
5773 ±60 K, and a pre-exponential factor, A, of 10 30±2 
molecules cm~ 2 s _1 . In §4, a comparison is made with 
previous results, from both thermodynamic and spec- 
troscopic techniques and the astrophysical implications 
of this work are discussed, focusing on applications as- 
sociated with dense molecular clouds. The conclusions 
are summarised in 55. 



2 EXPERIMENTAL METHOD 

The Nottingham Surface Astrophysics Experiment 
(NoSAE) was designed to measure surface binding en- 
ergies and sticking probabilities empirically and accu- 
rately, recreating the harsh interstellar environment un- 
der controlled laboratory conditions. A much fuller de- 
scription than can be given here is to be found else- 
where (Fraser, Collings & McCoustra 2001). Basically, 
the experiment consists of a stainless steel ultrahigh 
vacuum (UHV) chamber with an operating base pres- 
sure of 6xl0~ n Torr. The primary constituent (^90% 
) of this vacuum is H2, providing very similar con- 
ditions with the experiment chamber as those found 
in the ISM. The chamber is equipped with an effu- 
sive gas deposition system, quadrupole mass spectrom- 
eter (QMS), FTIR spectrometer and a Quartz Crystal 
Microbalancc (QCM). The QMS is used for tempera- 
ture programmed desorption (TPD) experiments (see 
§2.2), the FTIR for reflection-absorption infrared spec- 
troscopy (RAIRS), and the QCM for thin film mass 
determination. The grain mimic is an uncharacterised 
gold film surface, suspended at the end of a cryogenic 
cold finger, capable of reaching 7 K. The sample is 
radiatively heated from the reverse side by a halogen 
bulb. The sample temperature is measured using two 
E-type (Ni-Cr alloy vs. Cu-Ni alloy) and one KP-type 
(Au including 0.7% Fe vs. Chromel) thermocouples, po- 
sitioned on or near to the sample. The temperature can 
be controlled to better than 0.5 K, and measured ac- 
curately within 0.25 K. The upper temperature range 
of the system is 350 K. The whole system is linked via 
DAQ, GPIB, and RS232 communications to a pair of 
PC's, which are used to co-ordinate and control data ac- 
quisition from various instruments on the experiment. 
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2.1 Deposition 

H2O was obtained in situ from a liquid H2O sample that 
had been deionized, filtered for organic and inorganic 
solutes, and subjected to three freeze-thaw cycles under 
a vacuum of better than lxlO -7 Torr. H2O ice films 
were deposited at a rate of ca. 10 10 molecules s _1 on 
the gold film substrate held at a temperature of 10 K, 
from a quasi-effusive molecular beam of H2O, directed 
at 5° to the surface normal. At such low temperatures 
and slow deposition rates ballistic deposition results in 
the formation of high-density amorphous ice. Exposure 
of the substrate to gaseous H 2 was varied between 
1 and 100 Langmuir (L), (1 Langmuir is equivalent to 
lxlO -6 Torr s), and film mass was estimated using the 
QCM. An exposure of 20 L corresponded to a film mass 
of ca. 3.8xl0~ 7 g, which equates to a surface number 
density of H2O of 2.44xl0 16 molecules cm -2 , or a film 
thickness of ca. 0.03 /im, when proper account of the 
QCM sensor area is made. 



2.2 Temperature programmed desorption 

The sublimation characteristics of the H2O films 
were measured by temperature programmed desorption 
(TPD). In this method a linear temperature ramp is 
applied to the sample and the rate of desorption is mea- 
sured by monitoring the amount of adsorbate that des- 
orbs into the gas phase as a function of temperature. 
Desorption is an activated process which takes place at 
a rate given by 



dt s 



(1) 



where Agis the number of adsorbed molecules on the 
surface in molecules cm~ 2 , m is the order of the reaction 
and kd is the rate constant, given by 



Aexp 



Edes 

k B T 



(2) 



where A is the pre-exponential factor, Ed es is the bind- 
ing energy of a molecule on the ice surface in J, ks is 
the Boltzmann constant and T is the ice temperature 
in K. For convenience, in this paper the surface binding 
energies are expressed as Edes/^B, i.e. in K, and the 
units of A are determined from kd, i.e. for m — 1, A is 
expressed in s _1 , and where m = 0, A is expressed in 
molecules cm -2 s _1 . Equation (|l|) can be rewritten, to 
reflect the TPD signal that is actually measured during 
the experiment, i.e. 



dN s 
dT 



N m 



(3) 



where (3 is the heating rate, dT/dt. The recorded TPD 
signal peaks at some temperature maximum, T<j, which 
corresponds to the point at which the desorption rate 
from the surface is at a maximum, i.e. (d 2 Ns/dt 2 = 0). 



The TPD spectra are collected at increasing initial sur- 
face coverages, using the same linear heating ramp in 
each case. To a first approximation, the surface binding 
energy, Ed es , can then be calculated directly by substi- 
tuting for kdin equation (||) (using equation (0)), differ- 
entiating and equating to zero 



mC^cxp 



Edes 

kE/Td 



(4) 



Edes 

ksTj" 

provided the pre-exponential factor, A, is known. In a 
first order desorption process, where m = 1, A is as- 
sumed to be around 10 12 to 10 13 s _1 , approximately 
the vibrational frequency of the weak bond between 
the adsorbate molecule and the surface. It is also possi- 
ble to evaluate the order of the reaction from the TPD 
peak shape and peak maximum provided that the acti- 
vation energy for desorption and pre-exponential factor 
remain constant as a function of surface coverage. In 
practice however, desorption may not occur in a single 
step, the surface binding energy can vary across bind- 
ing sites on the surface or with surface coverage, and the 
pre-exponential factor, A, may differ by several orders 
of magnitude. A more rigorous result can be obtained 
by modelling the desorption system and using numerical 
fitting methods with a number of TPD spectra to eval- 
uate the values of A and Edes (see Woodruff & Delchar 
1986, Attard & Barnes 1998 for a more complete dis- 
cussion). 

In these experiments the ice layers were warmed at 
a rate of 0.02 K s _1 and the QMS, tuned to mass 18 
(H20 + ), was used to monitor the gas phase composition 
during desorption from the substrate. The QMS was 
operated in a line-of-sight configuration so that only 
molecules originating from the sample surface produced 
a signal at the detector (Fisher & Jones 1999). Sample 
coverages of 1, 2, 5, 10, 20, and 50 L (sub-monolayer to 
multilayer) were used. In all cases, no desorption was 
observed before 120 K and desorption was complete by 
170 K. 



3 RESULTS 

Fig. 1 shows a sequence of TPD spectra for H2O films 
deposited at 10 K, with increasing surface coverage. 
The spectra show the change in intensity of the mass 
18 (H 2 + ) signal at the QMS as the sample is slowly 
heated. An intense, asymmetric peak is observed at ca. 
140 K, associated with the desorption of the bulk ice 
film from the Au substrate. There is little or no change 
in the shape of the leading edge peak with increas- 
ing exposure. Furthermore, the peak temperature, Td, 
increases monotonically with exposure. These observa- 
tions are typical of TPD spectra in which the desorption 
process follows zero order kinetics, as illustrated by Fig. 
2. The figure shows results from a series of (a) zeroth 
and (b) first order TPD simulations (based on equation 
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Figure 1: TPD spectra for H2O ice films deposited at 
10 K on a Au substrate. Spectra are shown for surface 
exposures of 1, 2, 5, 10, 20, and 50 L. See text for details. 

3), at coverages of 2, 5, 10, 20, 50 and 100 monolay- 
ers of H2O, and under conditions similar to those in 
this experiment. In the zeroth order case (Fig. 2 (a)), 
the TPD curves are asymmetric, with coincident rising 
edges, and the peak temperature, Td, increases with ex- 
posure. This is typical of multilayer desorption, where 
the intensity of the TPD spectrum increases as more 
and more material is condensed onto the surface. The 
shift in Td (the peak position) occurs because the des- 
orption rate increases exponentially with temperature, 
so the rate can increase indefinitely until all the lay- 
ers have been stripped away, and the TPD signal falls 
rapidly to zero. The bonding between the layers resem- 
bles that in a condensed solid of the adsorbate, so the 
desorption energy, Ed es , and pre-cxponcntial factor, A, 
are the same for every layer. In contrast, the first order 
desorption curves are symmetric about a single T^ value, 
independent of coverage (see Fig. 2 (b)). This closely re- 
sembles monolayer desorption kinetics, when the nature 
of the bonding between the adsorbate and the underly- 
ing substrate is clearly important. In a multilayer sys- 
tem (e.g. icy mantles on dust grains, or the ice system 
studied here) the monolayer is generally more tightly 
bound to the underlying substrate than the subsequent 



adsorbate layers. Consequently, the monolayer desorp- 
tion peak appears at higher temperatures in the TPD 
spectrum than the multilayer desorption, reflecting a 
different Ed es value. Its line shape profile will also differ, 
reflecting both a different lvalue, and the different re- 
action kinetics. A more complete discussion of this com- 
plex relationship between TPD line shape profiles, bind- 
ing energies, and desorption kinetics can be found in 
Woodruff & Delchar (1986) or Attard & Barnes (1998). 
However, from a comparison between the experimental 
TPD spectra presented in Fig. 1 and the simulations in 
Fig. 2, it is reasonable to conclude that the results of 
this experiment are indicative of zeroth order kinetics. 
Furthermore, only a single peak is observed in the TPD 
spectra: even at low coverage it is not possible to decon- 
volve a monolayer desorption peak from the multilayer 
peak. Therefore no evidence exists to suggest that the 
desorption process is coverage dependant, nor that it 
is possible to distinguish between the multilayer and 
monolayer desorption processes. These observations are 
entirely consistent with previous reports on the TPD of 
ice from other hydrophobic, metal substrates (Kay et 
al. 1989; Dohnalek et al.1999, 2000). 

In some cases, and particularly at higher exposure, 
a slight shoulder is observed on the leading edge of the 
140 K TPD peak. Under these experimental conditions, 
the phase transition between amorphous and crystalline 
ice is known to occur between 120 and 140 K (Sack & 
Baragiola 1993). It might be expected therefore that 
all the desorbing ice should be in the crystalline phase 
by the time desorption occurs. However, given the rel- 
atively long timescale for this phase change (Dohnalek 
et al., 2000), it is likely that some desorption will occur 
from residual amorphous phase that has yet to undergo 
the phase change. The TPD peak of the amorphous 
ice will occur at a lower temperature than the main 
desorption peak since the vapour pressure of the amor- 
phous ice is around three times greater than that of 
crystalline ice (Kouchi 1987). This behaviour is clearly 
illustrated in Fig. 3. Traces (a) and (b) show the TPD 
spectra from H2O films of equivalent surface exposure, 
deposited at two temperatures, 10 and 100 K, producing 
high-density amorphous ice and low-density amorphous 
ice respectively. The films were cooled to 10 K and then 
treated as described in §2. Their TPD spectra are essen- 
tially equivalent, showing the amorphous ice desorption 
shoulder on the leading edge, which gradually flattens 
out as the crystallisation phase change dominates, and 
finally the major TPD peak, corresponding to desorp- 
tion of the crystalline phase. Two further films, with 
double the surface coverage of the amorphous films, 
were also investigated. Trace (c) shows results for an 
H2O film that was grown at 130 K then cooled to 10 K 
and warmed. The TPD spectrum from this ice film is 
narrower, with the onset of desorption occurring slightly 
later than the other spectra, and no evidence of the 
amorphous shoulder can be seen. These features are at- 
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Figure 2: TPD spectra from a series of (a) zeroth (A = 
10 30 molecules cm -2 s _1 ) and (b) first order (A— 10 13 
s _1 ) simulations (based on equation 3), at coverages of 
2, 5, 10, 20, 50 and 100 monolayers of H2O. In each 
case, the desorption energy, Ed es , was fixed at 5773 K, 
and the heating ramp, (3, was 0.01 Ks _1 . 



tributed to the purely crystalline nature of this ice film. 
Finally trace (d) in this spectrum, was produced by de- 
positing equal amounts of H2O at 100 and then 10 K. 
The total surface coverage in (d) is therefore equivalent 
to (c). The TPD trace for this sample shows evidence of 
the amorphous phase desorption, but the phase change 
appears to be much faster than in (a) or (b) , and spec- 
trum is dominated by desorption from the crystalline 
phase. Similar results have been found in previous stud- 
ies of ice desorption from various substrates, where the 
rate of crystallisation has also been determined (Speedy 
et al. 1996; Smith et al. 1996). 

A kinetic model was used to evaluate the surface 
binding energy, Ed es , and the pre-exponential factor, A, 
for the crystalline H2O ice system. To model the kinetics 
of the bulk ice desorption, a simple reaction scheme was 
constructed to describe the desorption process 



where the concentration of H 2 0(s) is equivalent to N$, 
the number of adsorbed molecules on the surface, and 
k p is the pumping coefficient of H2O in this system. A 
zero order rate equation, based on equation (]l|), was 
used to describe the desorption of molecules into the 
gas phase 

d[H 2 0(g)} 



dt 



fed 



(7) 



The rate coefficient, fed, is equivalent to that de- 
scribed in equation (0). A similar first order rate equa- 
tion was used to describe equation (||) 

d[H 2 0(g)} 



dt 



= -k d [H 2 0{g)] 



(8) 



The TPD spectrum then represents the temporal 
evolution of these coupled differential equations as the 
system temperature is raised in a linear manner. Us- 
ing an initial H2O ice surface density of 2.44xl0 16 
molecules cm~ 2 (corresponding to a 20 L exposure), 



■- 

cs 
- 

■- 

cs 



£ 



a. 
■- 
o 

a 

Q 

O 

(-1 

X 



(d) 



i'ii "Ih—imi " - V " % W * ■ »' ' I V< 1 1 1 "—" 



(c) 



(b) 



(a) 




J 1 I 1 I 1 L 



J 1 L 



H 2 0(s)^H 2 0(g) 



(5) 



20 40 60 80 100 120 140 160 180 
Temperature (K) 

Figure 3: TPD spectra for H2O ice films prepared (a) 
by 20 L exposure at 10 K, (b) by 20 L exposure at 100 
K, (c) 80 L exposure at 130 K, and (d) 20 L exposure 
at 100 K followed by 20 L exposure at 10 K. Each film 
was then cooled to 10 K prior to desorption. For clarity, 
spectra are offset and only shown from 100 K. 
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an initial temperature of 10 K, and a linear tempera- 
ture ramp of 0.02 K s _1 , model TPD spectra were cal- 
culated using a simple stochastic integration package *R 
(Houlc & Hinsberg 1995) and compared with the empir- 
ical data. The unknown parameters in the model, the 
surface binding energy, Ed es , and the pre-exponential 
for desorption, A, were then varied to best reproduce 
the empirical data. This was achieved for an A value of 
10 30±2 molecules cm -2 s _1 and an E des value of 5773 
±60 K. Fig. 4 shows the comparison of the model results 
with the empirical data. 
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Figure 4: A comparison between experimental and mod- 
elled TPD spectra. See text for discussion. 

It is much more complicated to evaluate the A and 
Edes values of the amorphous phase as the desorption 
and crystallisation processes are occurring simultane- 
ously. To evaluate the gas phase H2O concentration as 
measured in the experiment, it is necessary to couple 
equations describing the amorphous ice desorption rate 
and the crystallisation rate to equations ([j|) and (||). 
This work is currently ongoing, and will be the subject 
of a later paper. From the shapes and behaviour of the 



■f * Chemical Kinetics Simulator, Version 1.0, IBM, 
IBM Almaden Research Centre, 1995. Further infor- 
mation may hi obtained from Hie. CKS website at 

http://www.almaden.ibm.com/st/msim/ckspagc.html 



TPD curves in Fig.l and Fig. 3 however, it is clear that 
that the desorption mechanism for the amorphous ice is 
identical to that of the crystalline ice. Previous studies 
have drawn the same conclusion (Dohnalek et al., 1999, 
2000). Initial results from our model suggest that the 
^factor for amorphous ice is of the same order of mag- 
nitude, and that E ( i es is a few hundred K lower, than 
that of crystalline ice. 



4 ASTROPHYSICAL IMPLICATIONS 

The results obtained here compare favourably with pre- 
vious TPD measurements of H2O desorption under 
UHV conditions, from a range of surfaces, (Kay et al. 
1989; Haynes Tro & George 1992; Smith Huang & Kay 
1997; Chakarov & Kasemo 1998; Dohnalek et al. 1999, 
2000). When values for Edes and A are provided, the 
data are summarised in Table 1, where a comparison is 
also made with the spectroscopic results of Sandford & 
Allamandola (1988). With the exception of the latter re- 
sults, the table clearly shows that for both ordered and 
disordered hydrophobic surfaces, the desorption energy 
and rate constant are independent of the substrate. In 
each case, the desorption kinetics resemble zero order 
rather than first order reaction kinetics. 

None of the studies mentioned in Table 1 have iden- 
tified any differences between the desorption of H2O 
from H2O layers and the desorption of H2O from the 
underlying substrate. This indicates that H2O is only 
weakly bound to hydrophobic substrates, and that the 
adsorbate-substrate binding energy is comparable to 
the adsorbate-adsorbate hydrogen bonding. In astro- 
physical terms, this indicates that the same desorp- 
tion mechanism and desorption energy can be used to 
model H2O desorption from any hydrophobic surface, 
(e.g. silicon, graphite). Similar desorption studies on 
hydrophilic surfaces (such as metal oxides, quartz, and 
silicates) indicate that H2O bonds chemically to such 
surfaces. In such cases, the H^O-substrate system ex- 
hibits alternative desorption kinetics, with even higher 
activation barriers (for example Stirniman et al. 1996; 
Trakhtenberg et al. 1997; Hudson et al. 2001). However, 
these systems do exhibit identical H2O-H2O desorption 
processes in the regime where the ice surface is substrate 
independent ( i.e. the ice layer is thick enough that it 
is no longer influenced by the underlying substrate). 
Consequently, in ISM regions where relatively thick ice 
layers are formed, or the grain surface is hydrophobic, 
the underlying structure of the cosmic grain is not sig- 
nificant. Elsewhere, the grain material may influence 
the H2O desorption mechanism. It will not be possible 
to say conclusively which H2O desorption kinetics will 
dominate until the nature of interstellar grains has been 
clearly established. 

As Table 1 shows, the results obtained in this study 
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Table 1: A comparison surface binding energy, Edes, and pre-cxponcntial factor, A, measurements for (i) crystalline 
ice and (ii) amorphous ice. * units in s _1 not molecules cm" 2 s _1 as reaction was assumed to be first and not zero 
order, (a) This work, (b) Speedy et al. (1996), (c) Haynes et al. (1992), (d) Sandford & Allamandola (1988). 



differ significantly from those obtained by Sandford & 
Allamandola 1988. Firstly, the Ed es differs by ca. 800 
K for both the crystalline and amorphous ice phases, 
and secondly the desorption kinetics in this study are 
clearly zero order and not first order. These two parame- 
ters are both significant when modelling interstellar en- 
vironments, particularly when calculating the residence 
time of H2O on interstellar grains, or establishing the 
ratio of solid to gas phase material. 

Why do the two methods give such different results? 
The discrepancies must be related to the different as- 
sumptions that are made in interpreting the data under 
each set of experimental conditions. The most signif- 
icant difference between the two experiments is that 
the spectroscopic data are taken under isothermal con- 
ditions at base pressures of ca. 2xl0~ 8 Torr, whereas 
the TPD spectra are recorded in dynamic conditions 
at base pressures of 6xl0 -11 Torr. In the spectroscopic 
measurements, it was assumed that a first order rate law 
applied, and the value of .E^eswas determined indirectly, 
by monitoring the relative change in integrated intensity 
of a surface IR spectroscopic feature during isothermal 
desorption. Contrary to the assumptions made in Sand- 
ford & Allamandola's method, close to the desorption 



temperature, the sticking probability of H2O molecules 
on the surface is no longer unity, and over the time pe- 
riod of the experiment it is possible for re-adsorption 
to occur, both from vacuum contaminants and desorb- 
ing gas. The evaluation of isothermal measurements is 
also complicated by non-equilibrium effects, including 
both crystallisation and re-adsorption. In the TPD ex- 
periments it is possible to obtain information on the 
order of the reaction and Ed es directly, from the shape 
and position of the TPD curve. In UHV systems de- 
signed for TPD experiments, the pumping speed is so 
high that surface re-adsorption is not relevant. Unlike 
the spectroscopic data, the TPD spectra clearly show 
that crystallisation and desorption are occurring con- 
currently in the system. It is possible however to identify 
the amorphous, crystalline and 'phase change' regions 
of the TPD spectra from each other. With careful film 
preparation it is even possible to force the TPD spec- 
trum to exhibit only crystalline behaviour, as was seen 
in Fig. 3. Furthermore, Kay and co-workers have made 
a number studies using TPD methods to determine the 
rates of crystallisation and desorption when they are oc- 
curring at the same time, and have found no significant 
difference between the desorption rates they measure 
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Grain T (K) 



Half-life, t 1/2 , (yr.) 



w^r 



Zero order desorption kinetics, i.e. t\, — -^ 
where -/V S] o=1.15xl0 17 molecules cm~ 2 



First order desorption kinetics, i.e. t\, = ^ 



H2O on crystalline ice 
(a) (b) 



H2O on amorphous ice 
(a) (b) 



H2O on crystalline ice 
(c) (d) 



H2O on amorphous ice 
(c) (d) 
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6.8xl0 2m 
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IxlO 90 


2 xlO 102 


2.5xl0 84 


4.5xl0 106 


5.5xl0 9 ° 


1.3xl0 103 


1.6xl0 85 


30 


8.5 xlO 63 


2.1xl0 53 


3xl0 61 


3.5xl0 49 


4.5 xlO 64 


l.lxlO 54 


2 xlO 62 


2.2x10™ 


40 


8.4 xlO 42 


9.3xl0 34 


1.2xl0 41 


l.SxlO 32 


4.4xl0 43 


4.9xl0 35 


7.6xl0 41 


8.3xl0 32 


50 


2.1xl0 30 


9.1xl0 23 


6.7xl0 28 
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60 
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S.lxlO 15 


70 


8.4xl0 15 


24X10 11 


6.7xl0 14 


5.1xl0 9 


4.4xl0 16 
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4.3xl0 15 
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80 


2.7xlO n 


2.8xl0 7 


2.8xl0 10 


9.4xl0 5 


1.4xl0 12 


1.5x10 s 


1.8xlO n 


6.1xl0 6 


90 


8.4xl0 7 


24xlO i 


l.lxlO 7 


1.2xl0 3 


4.4xl0 8 


1.3xl0 5 


7.2xl0 7 


7.5xl0 3 


100 


1.3 xlO 5 


8.7X10 1 


2.1xl0 4 


5.6x10" 


7xl0 5 


4.6xl0 2 


1.4xl0 5 


3.6X10 1 


110 


6.8xl0 2 


8.7X10- 1 


1.3xl0 2 


7xl0~ 2 


3.6xl0 3 


4.6x10° 


8.1xl0 2 


4.5X1Q- 1 


120 


8.4x10° 


1.8xl0- 3 


1.8x10° 


1.9xl0~ 2 
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4.5xl0- 5 


4.2xl0" 2 


2.3xl0- 4 


1.4x10-2 


3.8xl0- 5 


150 


5.3xl0" 4 


4xl0~ 6 


1.5x10 4 


6xl0~ 7 
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2.1xl0~ 5 


9.4xl0~ 4 


3.8xl0~ 6 



Table 2: The half-life, ti/ 2 , of H2O molecules on H2O ice surfaces as a function of temperature under zero and first 
order kinetics. In each case, the half-life is defined in the table, and represents the time it takes for the surface 
population of H2O molecules on an interstellar grain to decrease to one half of its initial value, (a) Edes and A taken 
from this work, (b) Ed es taken from Sandford & Allamandola (1988), A taken from this work, (c) Edes taken from 
this work, A taken from Sandford & Allamandola (1988), (d) Ed es and A taken from Sandford & Allamandola (1988) 



isothcrmally or dynamically (Speedy et al. 1996; Smith 
et al. 1996; Dohnalck et al. 1999, 2000). Finally the 
experimental conditions under which the TPD spectra 
were obtained most closely resemble those in the inter- 
stellar medium. We therefore believe that these data 
represent the parameters for H2O that should be em- 
ployed in astronomical models. 

Table 2 shows the effect these new parameters have 
on the 'residence times' of H2O molecules on H2O ices 
and hydrophobic surfaces. It is not possible to define a 
residence time for a zero order process that is directly 
comparable (with the same units) with a first order 
process. Consequently, the results are presented here 
in terms of the half- life of the H2O surface population, 
rather than the 'residence time, 1/fcd' that is commonly 
used in the astronomical literature. The respective half- 
life values are defined in the table. It is important to 
note that the half-life of a zero order process also de- 
pends on the surface concentration of H2O molecules. 



In these calculations this value was fixed at 1.15 xlO 17 
molecules cm -2 , equivalent to around 100 monolayers of 
H2O ice, i.e. the number of H2O layers expected to ac- 
crete on a interstellar grain during the lifetime of a dense 
cloud (Hasegawa & Herbst 1993). The half- lives calcu- 
lated using the zero order ^4and Ed es values obtained in 
this work are given in column one (crystalline ice) and 
three (amorphous ice). The results are also compared 
to half-life calculations using the data from Sandford 
& Allamandola (1993). The half-lives calculated using 
Sandford & Allamandola's first order Aand Ed es val- 
ues are given in columns six (crystalline ice) and eight 
(amorphous ice). From Table 2 it is immediately clear 
that H2O can remain on the grain surface at signif- 
icantly higher temperatures than has previously been 
assumed, typically around 110 - 120 K rather than 90 - 
100 K. 
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5 CONCLUSION 

The surface binding energy and desorption kinetics of 
amorphous and crystalline H2O ice have been studied 
under conditions similar to those found in denser re- 
gions of interstellar clouds. These data have been used 
to calculate the residence times of H2O on H2O and 
hydrophobic surfaces as a function of temperature. The 
results imply that, in the dense interstellar medium, 
thermal desorption of H 2 ice will occur at signifi- 
cantly higher temperatures, and different rates, than 
previously had been assumed. These data are of funda- 
mental importance in the chemical modelling of many 
astrophysical environments. The effects may be partic- 
ularly pronounced in so-called hot cores, which are very 
dense clumps of gas, remnants of a collapsing cloud 
that formed a massive star. Irradiation of the clumps 
heats the cores to temperatures in the range 100 - 300 
K, when the ices evaporate, populating the gas with 
water and other trace molecules (cf. Millar 1993). The 
temperature of the evaporation may be crucial to relat- 
ing the observed chemistry of hot cores to the implied 
rate of warming of the central star (Viti and Williams 
1999). The results provide the astronomical community 
with reliable thermodynamic data pertaining to the des- 
orption of H2O under interstellar conditions. 
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